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The adiabatic manipulation of quantum states is a powerful technique that has 
opened up new directions in quantum engineering, enabling tests of fundamental con¬ 
cepts such as the Berry phase and its nonabelian generalization, the observation of 
topological transitions, and holds the promise of alternative models of quantum com¬ 
putation. Here we benchmark the stimulated Raman adiabatic passage process for 
circuit quantum electrodynamics, by using the first three levels of a transmon qubit. 
We demonstrate a population transfer efficiency above 80% between the ground state 
and the second excited state using two adiabatic Gaussian-shaped control microwave 
pulses coupled to the first and second transition. The advantage of this techniques is 
robustness against errors in the timing of the control pulses. By doing quantum to¬ 
mography at successive moments during the Raman pulses, we investigate the transfer 
of the population in time-domain. We also show that this protocol can be reversed by 
applying a third adiabatic pulse on the first transition. Furthermore, we demonstrate a 
hybrid adiabatic-nonadiabatic gate using a fast pulse followed by the adiabatic Raman 
sequence, and we study experimentally the case of a quasi-degenerate intermediate 
level. 


The precise control and manipulation of the states of 
multilevel quantum systems is a key requirement in quan¬ 
tum information processing. Using multilevel systems as 
the basic components of quantum processors instead of 
the standard two-level qubit brings along the beneht of 
an extended Hilbert space, thus ultimately reducing the 
number of circuit elements required to perform a com¬ 
putational task [T] . Superconducting qubits can be oper¬ 
ated as multilevel systems, and to date fundamental ef¬ 
fects related to three-level artificial atoms irradiated with 
continuous microwave fields - for example the Autler- 
Townes effect [H [3], coherent population trapping [3], 
and electromagnetically-induced transparency [5] - have 
been reported. An additional powerful concept is the 
adiabatic control of quantum systems, which in recent 
times led to new directions in quantum engineering, en¬ 
abling tests of fundamental concepts such as Berry phase 
[6] and its nonabelian generalization [7], the observation 
of topological transitions [5], and holds the promise of 
alternative models of quantum computation BM- 
In this paper we employ a superconducting circuit ir¬ 
radiated with two microwave fields to demonstrate a fun¬ 
damental quantum-mechanical process called stimulated 
Raman adiabatic passage (STIRAP) [T5l[T3]: the trans¬ 
fer of population between two states |0) and |2) via an 
intermediate state |1) using two control microwave fields, 
with zero occupation of the intermediate state at any 
time. The absence of population in the intermediate state 
is ensured by a two-photon destructive interference on 
the intermediate level which creates a dark state, while 
the adiabatic manipulation of the amplitudes of the con¬ 
trol fields rotates the dark state from the initial state 
|0) to the target state |2). The STIRAP technique was 
first demonstrated with sodium dimer molecular beams 
m and applied since then in many contexts in atomic 
physics. It has remained a topic of active research HH], 


and in recent years unexpected connections and analogies 
have been found. Recently there have been proposals on 
designing A systems using superconducting circuits with 
similar energy structure as atoms [siiin]. However, due 
to the present pre-eminence of the transmon as the build¬ 
ing block of the future superconducting quantum proces¬ 
sors, it would be most relevant to demonstrate this tech¬ 
nique on a typical device of this kind. The transmon |21] 
is a capacitively-shunted Cooper pair box strongly cou¬ 
pled to an electromagnetic transmission line resonator 
[22] , which allows the quantum non-demolition measure¬ 
ment of the state of the transmon by probing the res¬ 
onator’s resonance frequency with microwaves |23j . Here 
we aim at benchmarking the STIRAP protocol under re¬ 
alistic conditions - with a transmon that has finite deco¬ 
herence times associated to the excited states. 

In the transmon, it is not possible to transfer directly 
the population between the ground state |0) and the sec¬ 
ond excited state |2) due to the fact that the electric 
dipole momentum between these two states is vanishingly 
small. Thus this transfer must be done by first populat¬ 
ing the first excited state |1) using a tt pulse applied to 
the 0 — 1 transition, then transferring this population to 
the second excited state |2) using another tt pulse ap¬ 
plied to the 1 — 2 transition. This is called an intuitive 
sequence, and requires a precise control of the pulse se¬ 
quence, as well as of the frequencies of the microwaves. 
In contrast, when the transfer is realized using STIRAP, 
there is no population in the intermediate state during 
the process. Because of this, the decay associated with 
the intermediate state does not matter. The population 
transfer in STIRAP is based on adiabatically varying 
control microwave fields, applied in a counter-intuitive 
order, i.e. with the first pulse driving the 12-transition, 
and the second pulse, having a suitable overlap with the 
first one, driving the 01-transition. The adiabatic nature 
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of the population transfer makes the protocol very robust 
to the timing and the amplitudes of the pulses, as well as 
to the frequencies of the microwaves. Ideally, to ensure 
adiabaticity one should use long pulses. However, we find 
that this is not a strong restrictions, and the protocol can 
be run relatively fast using shorter pulses. 

We consider a transmon driven by two microwave 
pulses of amplitudes r2oi(t) and Hi 2 (f) with frequencies 
and (^^ 2 ^ (t), possibly slightly detuned from the 
corresponding qubit transition frequencies wqi and a;i 2 . 
In the dispersive regime and in the rotating wave approxi¬ 
mation with respect to the two driving tones (see Supple¬ 
mentary Information section I), the effective three-level 
Hamiltonian of the driven transmon reads 012 ! 


Hit) = 


n 
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0 f2oi(t) 0 

r2oi(t) 2^01 fli2(t) , 

0 f2i2(t) 2{Sqi + S 12 ) 


( 1 ) 


where the detunings ^oi = and 812 = W 21 — 

are defined with respect to the Lamb-shift renormalized 
transition frequencies of the transmon wgi = wgi -I- xoi 
and a;i 2 = W 12 -I-X 12 — Xoi; and for simplicity Hgi and Q .12 
are taken real. The rates xoi and xi 2 are the ac-Stark 
shifts of the corresponding transitions. If the two-photon 
resonant condition ^oi +<^12 = 0 is satisfied, the Hamilto¬ 
nian Eq. Q has a zero-eigenvalue eigenstate called dark 
state, |D) = cos0|O) — sin0|2) (see Methods), with 0 
defined by tan0 = Hoi(t)/f2i2(t). Due to the adiabatic 
theorem, by slowly tuning the drive amplitudes Hoi(t) 
and Hi 2 (t), such that tan0 goes from zero to infinity, 
the population transfer from the ground state to the sec¬ 
ond excited state can be realized. 

In the experiment we work with Gaussian pulses, 
parametrized as 




fll2(t) 


flgi exp 


O 12 exp 


’ 

2o-2 


( 2 ) 

(3) 


where ts is the time separation between the maxima of 
the pulses, a is the width of the pulse, and flgi: ^12 
are the amplitudes. We choose equal-amplitude pulses, 
rigi ~ Di 2 ~ H. The general condition for the adiabatic 
following is given by [TSJ[TS| |(±|(d/dt)|D)| ^ |tLi± — wdI. 
By using ([^ and performing a time integration (Supple¬ 
mentary Information section H), we obtain the global 
condition of adiabaticity for Gaussian pulses. 


> 1. (4) 

y/TT 


This global adiabatic condition shows that the transfer 
efficiency can be improved by making the pulses longer, 
which is in turn limited by the decoherence of the states 
|0) and 12). In principle, if a high transfer efficiency is de¬ 
sired, the shape of the pulses could be further optimized 

[21 US). 


EXPERIMENTAL RESULTS 

From the spectroscopy data as well as from the change 
of the Rabi frequency as a function of the detuning, 
we determine the Lamb-shift renormalized transition 
frequencies of the transmon ojoi/2tt = 5.27 GHz and 
oji 2 I 2 tt = 4.82 GHz. The dissipation rates extracted 
from independent measurements were Tig = 2.4 MHz, 
r 2 i = 5.2 MHz, and T^^ « T^’g = 0.4 MHz, while the 
Gaussian pulses were calibrated to Hgi/27r = 43.4 MHz, 
and ni2/27r = 38.2 MHz, with pulse width cr = 45 ns 
for both pulses (see Methods). This choice satisfies the 
adiabatic condition in Eq. Q by a factor of 4. 

Fig. [^) demonstrates the STIRAP protocol for our 
system, with the population of the state |2) reaching 
a maximum value during the pulse overlap and then 
slowly decreasing at large timescales due to the intrin¬ 
sic decay of the second excited state. The pulse se¬ 
quence is presented in Fig. The results from a 

quantum tomography measurement for three levels is 
shown schematically in Fig. &)• The traces are ob¬ 
tained by a homodyne measurement of the cavity re¬ 
sponse rj(T) = {{Ij{T)),{Qj{T))}, j = 0,1,2 for the 
reference states |0), |1), and |2), which are prepared by 
using fast microwave pulses applied to the 0-1 and 1-2 
transition (Supplementary Information section HI). Here 
Ij and Qj are the in-phase and the quadrature fields. 
For an arbitrary state p{t) at time t during the STI¬ 
RAP protocol the cavity response is a linear combina¬ 
tion [26] rmeas(t,r) = where pj(t) = 

Tr[p(t)|j)(j|] is the population of state |j). We extract 
the populations pj (t) using the Levenberg-Marquardt al¬ 
gorithm. We note also that the tomography method pre¬ 
sented here includes interrupting the STIRAP protocol 
before the full population is transferred, a protocol called 
fractional STIRAP Ezj. 

From the data we can see that the population in the 
second excited state reaches about 83% and the inter¬ 
mediate state gets slowly populated due to the decay 
12) —> |1). Thus, the main sources of nonideality are 
the decay and the dephasing of the upper levels, as well 
as imperfections in the timing of the short pulses used 
for calibration. A surprising finding is that the STIRAP 
transfer occurs in fact quite fast, predominantly in the 
overlap region: the optimal value is obtained not too far 
from the beginning of the wgi pulse. This means that 
the duration of the STIRAP process can be shortened 
by interrupting it earlier (with correspondingly shaped 
pulses) without too big penalties on the fidelity of state 
transfer, thereby reducing the overall gate length. 

Next, we present a complete experimental investiga¬ 
tion of this effect, together with numerical simulations. 
By varying the time separation between the two pulses, 
one can study the counterintuitive as well as the intuitive 
sequence. In Fig. 0 we present the occupation proba¬ 
bility of each of the states |0), |1), and |2) at different 
measurement times (horizontal axis), measured from the 
peak of the 01 pulse. The pulse separation (vertical axis) 
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FIG. 1. a) The state of the quantum system at different times 
during the adiabatic population transfer at fixed pulse sepa¬ 
ration. The dots represent experimental results, the dashed 
line is a simulation of the three-level dissipationless case, and 
the continues lines show the simulation with dissipation in¬ 
cluded. During the time evolution the intermediate state |1) 
remains almost unpopulated, b) Schematic of the pulse se¬ 
quence during the STIRAP. c) The population of the sys¬ 
tem at every time t is calculated from the cavity response 
{’’')} (only shown here), and the mea¬ 
surement data (cyan) is compared to the calibration data for 
the known states |0) (blue),|l) (green) and |2) (red). 


is the time interval between the peak of the 12 pulse 
and the peak of the 01 pulse. Thus, zero pulse sepa¬ 
ration means complete pulse overlap, negative values of 
the pulse separation refer to the counter-intuitive pulse 
sequence (a;i 2 is applied before woi pulse) whereas pos¬ 
itive values indicate the pulses are send in the intuitive 
order. The lower panels of Fig. show the correspond¬ 
ing simulation result using Lindblad’s master equation 
(see Eq. Q in Methods). From the experimental data 
we find that the optimal value for the pulse separation 
is « —90 ns, in close agreement with the result given 
by [24j and our own simulations. Importantly, one no¬ 
tices that for STIRAP the transfer efficiency is high and 
a relatively slowly varying function of the pulse overlap 
in a plateau extending from about -120 ns to about - 
80 ns, which demonstrates that accurate optimization of 
the pulse overlap is not required. In principle, popula¬ 
tion transfer can also be realized using the intuitive-pulse 
sequence, but in this case the method is sensitive to the 
timing of the pulses, as can be seen from the oscillations 
present at positive pulse separation in Fig. Clearly 
the stability plateaus of the STIRAP do not form in the 
case of the intuitive sequence. Moreover, in the presence 
of decoherence it is harder to achieve transfer fidelities 
as high as with counter-intuitive pulse sequence, because 
the pulse sequence cannot be reliably interrupted until it 
is completely finished. This robustness against errors in 
the timing of the pulses demonstrated here is one of the 
most significant benefits of STIRAP over non-adiabatic 
methods of population transfer. 

STIRAP is also very robust with respect to small er- 
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FIG. 2. Time-resolved two-pulse sequence as a function of the 
ujo\ and a;i 2 pulse separation, where 0 means that the pulses 
are completely overlapping. The upper figures show the ex¬ 
perimental result for the time evolution of states |0),|1) and 
12), whereas the lower row shows the corresponding simula¬ 
tion results, calculated using Lindblad’s master equation for 
a multilevel Jaynes-Gummings Hamiltonian (see Methods). 
The STIRAP sequence (blue rectangle) occurs at negative 
values of the pulse separation, while positive values of the 
pulse separation correspond to the intuitive sequence. 


rors in the frequencies of the drive signals. The effect is 
most notable along the axis where the two photon reso¬ 
nance condition is not violated, as demonstrated in Fig. 
1^ This makes the method resilient for example to fre¬ 
quency shifts that occur in qubits due to charge and flux 
fluctuations. Finally, STIRAP is not very sensitive to 
the widths of the pulses: we have tested cr = 40, 60, 70 
ns and obtained similarly large values for the transfer ef¬ 
ficiency. We hnd that when decreasing the pulse width 
below these values the adiabaticity condition is no longer 
fulfilled and oscillation start to appear, while if it is too 
large the effect of decay becomes significant. 

During the operation of a quantum processor, it might 
be efficient to use both fast pulses and adiabatic pulses. 
Thus it is interesting to look at hybrid pulses, where for 
example a fast pulse is applied to the transition 01 before 
a STIRAP sequence, see Fig. This corresponds to 
starting the adiabatic population transfer protocol from 
a superposition between the ground state and the first 
excited state. During the evolution of the system, this 
results in a superposition between the dark state and the 
intermediate state, as shown in Fig. Towards the end 
of the sequence, the population on level 2 is stabilized 
(to a reduced values) but the populations of states 0 and 
1 oscillate due to the Rabi coupling provided by the last 
woi pulse. 

In Fig. we show how the STIRAP can be re¬ 
versed by applying one additional Gaussian pulse at the 
Wi 2 transition. The reversal of the STIRAP is essen¬ 
tial in creating adiabatic gates, for example to realizing 
a single-qubit phase gate [3S]. The idea is to hrst ex¬ 
cite the system to state |2) using STIRAP, and then di- 
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FIG. 3. STIRAP as a function of the detunings of the drive 
pulses. The offset from zero of the sum of the detunings (5oi + 
5i2 corresponds to the violation of the two photon resonance 
condition, to which the transfer is mostly sensitive. Nonzero 
values of the difference of the detunings 5oi — 5i2 correspond 
to the violation of the single photon detuning, against which 
the system is robust. The lower left corner is the measured 
ground state population, while the corresponding simulation 
results are shown in the other plots. 




a) 

„ 200 
1 150 
glOO 

i 50 
? 0 
§200 
“ 150 

CO 

■5 100 
CL 

50 


0 200 -200 0 200 -200 0 200 400 

Time [ns] 


1 


o 

13 


a. 

o 

Q. 


0 


-200 -100 0 100 200 300 400 


Time [ns] 


FIG. 5. STIRAP reversal: three adiabatic pulses are used to 
first excite the system to the second excited state and then 
bring it back to the ground state, a) The time evolution of the 
state (with the blue rectangle indicating the region where the 
process is the most efficient). The upper row is experiment, 
while the lower row is simulation, b) Schematic of the applied 
pulse sequence. 
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FIG. 4. a) The STIRAP efficiency depends on the initial 
population of the ground state. The Rabi pulse before the 
STIRAP pulse sequence transfers some of the population to 
the intermediate state |1), which reduces the population of 
level 12) at the end of the pulse sequence. The upper row 
shows the measurement results for the three states. The cor¬ 
responding simulation results are shown on the lower row. b) 
The pulse sequence used to generate the measurement results. 


rectly reapply STIRAP to bring the system back to the 
ground state. This demonstrates that STIRAP can also 
be used to transfer population from state |2) back to 
ground state. 

Finally, an interesting situation is that of a degenerate 
or quasi-degenerate intermediate state. To realize this 
degeneracy we did several thermal cycles of the sample, 
until we created a defect in the insulating layer of the 
junctions, observed by the splitting of the spectral line 
into states |la) and |lb). The formation of this type of 
defect has been observed recently in flux qubits as well 
P5] . In this case an interesting two-blob structure can be 
seen in the detuning plot, as shown in Fig. where each 
blob corresponds to the resonance conditions of either 
|la) or |lb). 


DISCUSSION AND CONCLUSIONS 

We have demonstrated the adiabatic transfer of pop¬ 
ulation from the ground state |0) to the second excited 
state 12) in a three level transmon circuit. This bench¬ 
marks the STIRAP protocol as a valid procedure for 
quantum information processing applications in circuit 
QED. As an added bonus, we note that while in the pre¬ 
vious continuous-wave experiments on various types of 
qubits [2H5] it was not possible to asses unambiguosly 
the effect of quantum interference from the spectrocopic 
data Uni, our experiment establishes it conclusively. We 
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FIG. 6. The splitting of the first excited state to the states 
I la) and |1(,) suppresses the transfer efficiency when the one 
photon detuning is between the split levels. The energy sep¬ 
aration between the levels is 15 MHz. 

also studied the effect of having a finite initial state pop¬ 
ulation in the dark state, we demonstrated the adiabatic 
reversal of the STIRAP protocol, and we realized the 
transfer of population also via a quasi-degenerate inter¬ 
mediate state. The main limiting factor is decoherence, 
which however for superconducting qubits is improving 
fast with new designs and fabrication techniques. We 
find that the process is very resilient to errors in the tim¬ 
ing of the pulses and the frequencies of the microwave 
drives. 


METHODS 


For zero-detuning Jqi = 0 we have $ = 7r/4. The adi¬ 
abatic tuning of the pulse amplitudes ensures that the 
state of the system follows the |D) eigenstate instead of 
being nonadiabatically excited to either one of the |±) 
eigenstates, which contain contributions from the state 
| 1 )- 

The transition frequencies of the transmon (I’oi/27r = 
5.27 GHz, £l’i2/27r = 4.82 GHz are obtained from the 
spectroscopy data as well as from the Chevron patterns 
of the detuned Rabi oscillations. The measurements of 
Rabi oscillations on the 1-2 transition are done by first 
applying a tt pulse to the 0-1 transition followed by the 
Rabi pulse on the 1-2 transition. To obtain the decay 
rate Tlo we used a tt pulse on the 0—1 transition and 
measure the population of the level 1 as it decays. To 
find ri 2 we use a tt pulse applied to the 0 — 1 transi¬ 
tion, then a tt pulse on the 1 — 2 transition, and we fit 
the data with a three-level model which includes the de¬ 
cay 2 —>■ 1 and 1 —?► 0. We obtain the energy relaxation 
rates Tig = 2.4 MHz and r 2 i = 5.2 MHz. From Ramsey 
interference experiments we get F^^ « F^q = 0.4 MHz. 
The Gaussian pulses are created by a high-sample rate 
arbitrary waveform generator followed by mixing with a 
microwave pulse (see Supplementary information). They 
are calibrated using the Rabi data, resulting in Hoi/27r = 
43.4 MHz, and fli2/27r = 38.2 MHz, with standard de¬ 
viation cr = 45 ns for both pulses. The quantum to¬ 
mography procedure follows Ref. (see Fig. and 
the Supplementary information for details). In our ex¬ 
periment, the tomography pulse interrupts the STIRAP 
protocol before the full population is transferred, thus it 
could be extended as well to the investigation of arbi¬ 
trary superpositions between the ground state and the 
second excited state in the so-called fractional STIRAP 
protocols [27]. For the simulations we use the three-level 
Lindblad master equation 


The effective Hamiltonian of the system in a doubly 
rotating frame and under the rotating-wave approxima¬ 
tion takes the form Eq. Q. If the two-photon de¬ 
tuning (5oi + Si 2 is zero, the ei genvalues of the above 
Hamiltonian are a;_|_ — <^01 + \/— 
Sqi — + -I- 11^2 j wd = 0, with correspond¬ 

ing eigenvectors 

|-|-) = sin$|B) -I- cos$|l), 

|—) = cos<i)|B) — sin$|l), (5) 

|D) = cos0|O) — sin0|2), 

where the zero-eigenvalue is called the dark state, and 
the state orthogonal to it |B) = sin0|O) -I- cos0|2) 
is called the bright state. The angle 0 is defined by 
tan0 = Hoi(t)/Hi 2 (t) and it parametrizes the rotation 
in the {|0), |2)} subspace (spanned also by the bright and 
dark states), while 


dt^ h 


p,H 


T -^rel [p] T ■^deph [p]: 


where the Lindblad superoperators are 


•^rel[p] — (2(Toip(Toi — (Tup — pUll) 


+ ^- (2cri2pcr2i — cr22P — pcn) 


for the relaxation and 

-Cdephip] = ( 2 (ToiP(Toi - CTiip - pern) 

pV? 

+ ^ (2cri2pcr2i - cr22P “ 


( 7 ) 


( 8 ) 


( 9 ) 


$ = tan ^ 


_ Vi^oi + ^ 12)72 _ 

V (^01 + ^12 + '^0 l )/2 + ^ 01 / 


(6) for the pure dephasing term and <7^ = |i)(j|. 
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Supplementary material 

I. SAMPLE AND MEASUREMENT SETUP 

The superconducting QED circuit [21] consists of a 
capacitively-shunted Cooper pair box (a transmon), with 
tunable energy level separation through the application 
of an external magnetic field. The qubit can be excited by 
applying a continuous wave or pulsed microwaves. The 
readout system is realized as a A/4 coplanar waveguide 
resonator which allows the quantum non-demolition mea¬ 
surement of the state of the transmon [221 [23] . A picture 
of the sample is shown in Fig. Izl^) ,b), and a schematic of 
the measurement setup at room temperature and in the 
dilution refrigerator is presented in Fig. &)• 

The transmon is essentially an artificial atom with en¬ 
ergy levels hiOj yielding transition frequencies denoted 
by Wjj+i = iOj+i — ujj, and can be modelled as an an- 
harmonic oscillator. With the standard notations in the 
field of superconducting devices, Ej denoting the Joseph- 
son energy of the qubit and Eq the total charge en¬ 
ergy (including the shunt capacitor), the anharmonicity 
is fiLJi 2 — fiujQi « —Eq in the asymptotic limit Ej ^ Ec- 
The usual way to manipulate the state of the circuit QFD 
system is to apply microwave drive signals either to the 
gate of the transmon or to the coupled transmission line 
cavity. If the frequency of the drive matches the 

transition frequency between two states of the system, 
the system goes through Rabi oscillations with frequency 
fljj+i, resulting in transfer of population between the 
two states. 

We model our system using the driven many-level 
Jaynes-Cummings Hamiltonian coupled to a resonator of 
frequency Wres ■ A measurement tone of frequency Wmeas is 
coupled with strength Cmeas into the resonator from the 
input capacitor of the coplanar waveguide cavity. The 
total Hamiltonian is then 

N 

Hjc = ^ Hujj\j){j\ + hwresa'^a 

3=0 

N-1 

+ Ho+I (Ij + 1)01® + h-c.) 

j=0 

3=0 

+ ^Cmeas (aie-““““* -f h.c), 

where |j) are the eigenstates of the Wlevel transmon 
corresponding to the eigenenergy hujj. Here we have ap¬ 
plied the rotating wave approximation to the driving and 
the measurement fields, retaining only energy-conserving 
terms. We also only keep the coupling between consec¬ 
utive energy levels with the corresponding near-resonant 
fields. 

The single mode resonator, described by the annihi¬ 
lation (creation) operator a (d’^), couples to the qubit 




FIG. 7. a) Optical microscope image (false colors) of the main 
parts of the sample, b) SEM image of one of the Josephson 
junctions in the transmon. c) Simplified schematic of the 
room temperature electronics, and the cryogenic setup and 
wiring used in the experiment. 


transition j j + 1 with the Jaynes-Cummings coupling 
strength g^j+i- 


In the dispersive regime, when the qubit is de¬ 
tuned from the cavity and the number of pho¬ 
tons in the resonator is not too large such that 
4(d'^d) — (jjj — Wres)^] "C 1, we can perform 

a generalized Schrieffer-Wolff transformation with dis¬ 


placement operator D = exp 


A 


where 


N 


^ = E 


ffjj+i 


j^Q OJj LOj-es 


(d|j +1)01-h.c.), (11) 


on the Hamiltonian in Eq. (10), resulting in El 
DHD'f. Next, we use 


e^He-^ = H 


A,H 


+ I [^’ L- 


A.H 


( 12 ) 


with gjj+i/(ujj^i—ujj—ujres) as a small parameter and re¬ 
taining only the first-order terms. This results in the can¬ 
cellation of the Jaynes-Cummings interaction between 





















































































the qubit and the resonator, 

N 

i?disp-JC ^ + Xj-I,j)\j){j\ 

i=i 

— ^Xo,id'*'d|0)(0| + hujj-es^^^ 

N 

+ (Xi-I J - Xj,j+i) a^a\j){j\ 


i=i 

N-1 




i=o 


(S2) 


j + l)(j|e + h.c. 


+ femeas (o^e + h.c) . 


Here the dispersive shifts Xjj+i s-re defined as 
Xj,j+i = 


n2 

9j,J + 1 


UJj 


(13) 


(14) 


which transforms the Hamiltonian as H —>■ UHIJ^ + 
iH{dU/dt)W. We also perform a rotation with respect 
to the cavity at the measurement frequency, C/meas = 
exp [zWineasQ^a^] resulting in a similar transformation 
H + ih{dU^eB^s/dt)U^^^^. Finally, this 

yields 


^disp-JC — X! ^ ( ^3 ^0 + Xj-i.i X! ^fc.fc + l j \d){j\ 

j=l \ k=0 / 

- %O,ia'''a|0)(0| +h{uJres - Wmeas) d'I'd 
N 

+ XI ^ (Xj-Ij - Xjj+i) a^a\j){j\ 


i=i 

AT-l 
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(|j + 1)01 +h-C.) +h€mea.s {d + O^) ■ 


3=0 


In obtaining Eq. (131, several terms have been neglected. 


A two-photon Hamiltonian resulting from the Jaynes- 
Cummings interaction couples states separated by two 
ladder indices. 


{2ujj+i - UJj - ujj+2)dgjj+igj+ij+2 -2|^-_^2)(j|+h 


N-2 

V . 

■ _Q ‘^{ojjj.2 OJjj-i 6di-es)(^j-t-1 OJj Wres) 

(15) 

Then there are terms resulting from the commutator of 
A with the drive. 


E 


93,3 + 1 


UJj+i 


— UJi — CJr, 


d -i, +, 

^y^e-*“-^-‘|j + l)0-l| 


- + 2){j\ + 

(|j + l)(j + l|-|j)0l)], 

and with the measurement pulse 


(16) 


(19) 

Exciting the qubit on the energy level j > 1 results in 
a measurable shift of the resonator transition frequency 
by Xi-i,i + Xo,i ~ Xj,j+i^ thus enabling a quantum non¬ 
demolition measurement of the transmon state. This is 
realized by monitoring the response of the cavity near the 
cr.esonance under the application of a measurement pulse. 
Before the measurement held is applied, that is during 
the STIRAP pulses, the number of photons in the cav¬ 
ity, and therefore the ac-Stark shift on the transmon is 
negligible. In this case the only effect of the resonator 
on the qubit is that its vacuum huctuations produce a 
small Lamb-shifted renormalization of the energy levels, 
UJj = UJj + Xi-i,i foi' J ^ 1 uud cuo = ujq, with corre¬ 
sponding transition frequencies cujj+i = ujj+i — ujj. Let 
us introduce the detunings as 


A Xjj'-i-i ojj X3-i,3 


he„ 


93,3 + 1 




[|j + l)(j|e -I- h.c.] . 

(17) 

Several arguments can be invoked to neglect these terms: 
some are second-order detuned processes such as Eq. 
(15), and the rest would produce only hrst-order cor¬ 


rections to the dominant driving term. Moreover, when 
the qubit is driven, the resonator is not populated, and 


therefore the relevant terms left in Eq. (16) contribute 


only as a renormalization of the decoherence of the sec¬ 
ond excited state into the resonator. 

Next, we set the ground state as a reference for the 
energy levels by subtracting a quantity X)Xi |j)(j | 
from the resulting Hamiltonian. Eurther, we move into a 
multiple-rotating frame defined by the transformation 


N 


m = E 


exp 




Wfc,fe+ii 






for j > 1 and i5oi = uii + xoi — wq — uj^'’ for the first 
transition. 

Next we focus on the STIRAP by considering only the 
three lowest eigenstates of the transmon, i.e. the states 
|0), |1), and 12). This eventually leads to a very simple 
three-level form for the system Hamiltonian 


H(t) 


0 Hoi(t) 0 

iloi(t) 26oi ^i2{t) 

0 Hi2(t) 2((5oi-I-(ii2) 


( 20 ) 


To simplify the notation, in the expression above as well 
as in the main paper we have eliminated the comma be¬ 
tween subscript indices whenever this does not lead to 
confusion, for example we write 50,1 = 9oi, <^ 0,1 = 1^01 j 
^^0,1 = f^oi- 


(18) 
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II. DERIVATION OF THE GLOBAL 
ADIABATICITY CONDITION 


Adiabatic processes should satisfy |(±|((i/dt)|D)| «C 
|a;± — wdI, which, when employing ([^ leads to 






< 1, 


( 21 ) 


for (5oi = 0 and sin$ = 1. 

For Gaussian pulse shapes with equal amplitude 

+2 


Ooi(0 = ^ 

0 i 2 (t) = exp 


t- 

2a2 


( 22 ) 

(23) 


we can find a global adiabatic condition by integrating 
Eq. (21) over time 


OO 

— ^ y A/Ooi(t)^ + Vti2{tY dt. 


(24) 


The remaining integral does not have an analytical solu¬ 
tion. However, it can be shown that it is an increasing 
function of moreover, the values at tg = 0 and tg = oo 
can be calculated analytically, and they are 2y^crO and 
2-\/^crH, respectively. Thus, we have to select the mini¬ 
mum value 2-y/7r(TH to obtain the constraint, and we get 

4 

-^aVL > 1. (25) 


III. CALIBRATION AND MEASUREMENT 
PROTOCOL 

We first determine the transition frequencies of the 
transmon by measuring the Rabi frequency of a driven 
transition as a function of the drive frequency. This gives 
woi/Stt = 5.27 GHz, Wi2/27r = 4.82 GHz. 

The used measurement setup is shown in Fig. To 
characterize the STIRAP process we determine the state 
of the three level transmon by probing the coupled trans¬ 
mission line cavity with a measurement signal. Due to 
the dispersive shift in the cavity resonance frequency, 
the response of the cavity to the probing microwave sig¬ 
nal is different depending on the state of the transmon. 
To measure the cavity response we apply homodyne de¬ 
tection scheme, where the incoming microwave signal is 
downconverted to DG using an IQ-mixer. The resulting 
in-phase and quadrature signals are given by 

{i(j)) =-V{a + a^) 

( 0 ( t )) =ir]{a-d^), 


the corresponding cavity response, we can associate the 
response of each of the states with a variable rj (r) = 
(''“))) {Qji^))} then describe the cavity response 
of every other state as a linear combination of Vj (t) 

?’meas(t,'r) = ^ Pj{t)rj{T). (27) 

i=0,1.2 

This is a linear combination of the responses of each state, 
weighted by the population Pj{t) — Tr[p{t)\j){j\] of the 
state, J2j=o 1 2 Pji^) ~ 1 — Pji^) — 1 

be understood by introducing the operator corresponding 

to r^easit^r), fi„eas(r) = Z)i=o.i,2 G ('^)Ij) 01 ^nd thus 

obtaining Eq. (|^ as r^eas{t,T) = Tr[p(t)fmeas(T)]. 

For a given trace in the I — Q plane, we 

determine poi Pi, and p 2 by applying the least square 
fit method to invert Eq. ( p^ . To implement the least 
squares fit method we employ the Levenberg-Marquardt 
algorithm, using the measurement response data up to 
r = 1500 ns. As pre-processing, the data is multiplied 
with an exponential weighting function exp[—r/w], with 
w = 700 ns. This allows us to use predominantly the 
data from the beginning of the response curves, where 
the effect of decoherence on the distinguishability of the 
states is minimal. 

In Fig. Ic) we present the (/(r)) traces correspond¬ 
ing to the system in the ground state (blue), first excited 
state (green), and second excited state (red), together 
with the trace for the density matrix at a time t = 450 
ns (cyan) during the STIRAP protocol. The calibration 
trace of the first excited state is obtained by applying 
a resonant tt pulse to the |0) —?► |1) transition, while 
the calibration trace of the second excited state is ob¬ 
tained by first populating the first excited state with a 
TT pulse, and then applying another tt pulse in resonance 
with the |1) —>■ 12) transition. This allows us to infer the 
density matrix of the system, given a set of measured 
, {Q{t)) traces. The main source of error in this 
procedure is produced by decoherence. The decay from 
the upper states during the calibration pulses results in a 
small inaccuracy, which produces a small overestimation 
of the populations in the tomography. We minimize this 
error as postprocessing by including the effect of deco¬ 
herence in the calibration. 


IV. STIRAP AS A ROUTE TO HOLONOMIC 
QUANTUM COMPUTING 

The phases of the driving fields, which we ignored so 
far, can lead to interesting new ways of qubit manipu¬ 
lation through the accumulation of geometric phases. If 
the phase factors are retained in the driving fields, we get 
the following three-level form for the system Hamiltonian 


where p is a. factor describing the losses of the conver¬ 
sion and the other constants. By preparing the trans¬ 
mon in the states |0), |1), and |2) and then measuring 


Hit) 


h 

2 


0 

Qoi 

0 


“01 

2(5oi 

0*2 


0 

Oi2 

2(<5oi -I- S 12 ) 


(28) 
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Let us denote by ip the phase difference between the driv¬ 
ing fields, therefore we can take fioi = |f^oi|, ^12 = 
|Ooi| exp(—At two-photon resonance i5oi -I- 612 = 0, 
the eigenvalues/eigenvectors problem at a time t has 
the solution u}+ = <5oi -I- \/Sqi + |floiP + |^^i 2 p, = 

(5oi - -y/cJoi + + |f^i 2 p, and wd = 0, with corre¬ 

sponding eigenvectors 


nection formula 


fRj 

7Berry = W (D|VR|D)(iR, (30) 

Jr, 


where the integral is taken along a contour in the (0, ip) 
space, 


R = 



(31) 


|-|-) = sin<i)|i?) -I- cos$|l), 

|-) = cos$|B) - sin$|l), (29) 

|D) = cos0|O) — sin0e*'^|2). 


where the bright state is |i?) = sin0|O) -I- cos 06*“^ 12). 
The angle 0 is defined by tan0 = |floi(i)|/|f^i 2 (i)| 
and parameterizing the rotation in the {|0),|2)} sub¬ 
space, and $ is defined in the same way, as an angle of 
a right triangle with vertices (IfloiP + |f^i 2 p )/2 and 

\/(|f^OlP + |^^12|2+<5o\)/2 + S 01 /V 2 . 


For a general trajectory starting at an initial time ti, this 
yields at time t [2H] 

rvit) 

7Berry(i) = - / Sm^ Qdip. (32) 

Jviti) 

Because the energy eigenvalue of the dark state is zero, 
from the time ti to time t the system would only pick up 
a geometrical phase, 

|D(t,)) ^ e*T'®-">'(*)|D(t)). (33) 

After a full STIRAP process the state then changes as 

|0) -J. _e0Be„y(t/)+»v(*/)|2). (34) 


The Berry phase accumulated during a STIRAP se¬ 
quence can be calculated from the standard Berry’s con- 


In holonomic quantum computing, the geometric phase 
7Berry IS uscd to construct single-qubit phase gates. 






